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ABSTRACT 

Existing studies of the atomic hydrogen gas content in distant galaxies, through the absorp- 
tion of the 21 -cm line, often infer that the total column density, A/hi, is anti-correlated with 
the linear extent of the background radio source, d cm . We investigate this interpretation, by 
dissecting the various parameters from which Nm is derived, and find that the relationship 
is driven primarily by the observed optical depth, r b s , which, for a given absorber size, is 
anti-correlated with d cm . Therefore, the inferred iVni — d em anti-correlation is merely the 
consequence of geometry, in conjunction with the assumption of a common spin tempera- 
ture/covering factor ratio for each member of the sample, an assumption for which there is 
scant observational justification. While geometry can explain the observed correlation, many 
radio sources comprise two radio lobes and so we model the projected area of a two com- 
ponent emitter intercepted by a foreground absorber. From this, the observed r bs — 4m 
relationship is best reproduced through models which approximate either of the two Fanaroff 
& Riley classifications, although the observed scatter in the sample cannot be duplicated us- 
ing a single deprojected radio source size. Furthermore, the trend is best reproduced using an 
absorber of diameter ~ 100 — 1000 pc, which is also the range of values of d mn at which the 
21 -cm detection rate peaks. This may indicate that this is the characteristic linear size of the 
absorbing gas structure. 
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1 INTRODUCTION 

The 21 -cm transition of hydrogen (Hi) probes the cool, star- 
forming, component of the neutral gas throughout the Universe. 
Due to the low probability of the transition, compounded by the in- 
verse square law, this is difficult to detect in emission at redshifts 
of z £ 0.1 with current instruments, although in absorption the line 
strength depends only upon the total hydrogen column density and 
the strength of the background continuum (see Eq. [T}. At z > 0.1, 
this transition has been detected along eighty different sight-lines, 
half of which occur in galaxies that inter vene the sight-lines to 
more distant quasars (see lGupta et al.l2012l and references therein), 
with the other half occuring within the host g alaxy of the radio 
sourc e itself (associated absorption; compiled in Curran & Whiting 
20101 with more recent results published in lCurranetal.ll201lj |bl 
2013al lbh. 



Of the associated absorbers, detection rates appear to be 
higher in those objects classified as "compact"Q Furthermore, 



E-mail: sjc@physics.usyd.edu.au 

That is, gigahertz peaked spectrum sources (GPS, ^ 



IPihlstrom et alj d2003b reported an anti-correlation between the de- 
rived atomic hydrogen column density, JVhi, and the projected lin- 
ear size of the background radio emission, d Gm , a finding which 
has since been sustaine d t hrough the addition of further data (e.g . 
Gupta & Saikial l2006allbt iGupta et al.1 120061 : lOrienti et al.1 l200d: 
Chandola et al.ll201 lh . 



This interpretation suggests that the smaller the radio source, 
the denser the absorbing medium, which could be consist ent with 
either the "frustration scenario" dvan Breugel et al.lll984l). wher e 



1 kpc) and com- 



the jets are confined, or the "youth scenario" dFanti et al.lll995l), 
where the gas has yet to be expelled (supported by IPihlstrom et al] 
120031) . In either case, there is one major flaw in the interpretation 
that the column density is anti-correlated with the jet size: It is the 
observed velocity integrated optical depth of the 21 -cm line that 
is actually measured, rather than the H I column density (Fig. [TJ 
top), which is derived by assuming a spin temperature and cover- 



pact steep spectrum sources (CSS, ~ 1 — 10 kpc), in addition to their 
various sub-classes - compact symmetric objects (CSO), compact flat spec- 
trum (CFS), high frequen cy peaker galaxies (HFPs) and low power com pact 
(LPC) radio sources (e. dO Deali998 ; Fanti 20()o HOrienti et al.l2O06l) . 
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Figure 1. The velocity integrated optical depth (top), the peak observed 
optical depth (middle) and the full-width half m axima (bottom) of the H I 
21-cm line versus the projected linear size. From lPihlstrom et alj fe003l) . 



ing factor (usually T sp i n = 100 K and / = 1). Where column 
densities are knowrQ the observable quantit y T avln _/ f is seen to 
vary by a factor of at least ~ 170 (fro m 60 K. ICurran et alj|2007l to 
9950 K. lKanekar & Chengaluj[200l) . and so this assumption can- 
not be justified. Furthermore, due to the scatter i n the full- width 
half maxima of the 21-cm profiles (Fig.[TJ bottoml. lPihlstrom et alj 
acknowledge that the correlation is driven by the observed optical 



Usually in z a b s ~ 1.7 damped Lyman-a absorption systems (DLAs). 



depth (Fig. [TJ middle). This therefore suggests that the observed 
relationship arises directly from an optical depth-linear size anti- 
correlation, the possible cause of which we explore in this paper. 



2 OPTICAL DEPTH AND RADIO SOURCE SIZE 

2.1 Observed optical depth and covering factor 

The atomic hydrogen column density along the line-of-sight, Nhi, 
is related to the velo city integrated optical de pth of the 21-cm ab- 
sorption, Jrdv, via dWolfe & Burbidgell 19751) : 



N m = 1.823 x 10 i!i T s , 



r dv , 



(1) 



where T sp i n [K] is the mean harmonic spin temperature of the gas 
and the optical depth is given by 



t = - In 1 - 



7"obs 

T 



(2) 



Here r D b s = AS/S is the observed optical depth of the line, 
where AS is the spectral line depth and S the continuum flux 
of the background source. The covering factor, /, quantifies how 
effectively the absorber intercepts the flux (ranging from zero to 
unityEl meaning that r bs ^ T. In the optically thin regime, 
Eq. [2] simplifies to r ~ T b s / f, so that Eq. [JJis approximated by 
N m « 1-823 x 10 18 (Tspin//) fr obs dv. Therefore, with knowl- 
edge of the value of T sp i n / f, the column density can be obtained 
from the observed integrated optical depth. However, this ratio is 
generally unknown^, with the assumption of a sin gle value being 



applied to yield the column d en sity in each source ( Pihlstro m et al 



2011). 



20031 : iGupta & Saikiall2006allbl iGupta et alj|20o3 



Chandola et al 



We now consider, without evidence to the contrary, that the 
column density of a foreground absorber is independent of the frac- 
tion of background radio flux it covers. In the optically thin regime, 
for a given absorber the observed optical depth is proportional to 
the covering factor (Eq.|2}, which in turn is anti-correlated with the 
size of the background source. That is, r b s oc / = (A^hs/Aem) , 
where the last term is the ratio of the cross-sectional area of the 
absorber to th at of the background emitter (where z a b s ~ z C m, cf. 
ICurranll2012l) . 

Therefore, based upon its definition alone, the observed opti- 
cal depth is anti-correlated with the size of the radio source. This is 
seen in a plot of these two quantities (Fig.[JJ middle), with the scat- 
ter being due to differences in the actual optical depths, r, and, 
possibly, the absorber sizes. In Fig. [2] we re-plot this and over- 
lay the variation of r bs with d cm for r b s = ,f{l — e" T ) and 
/ = (^absMem) = (dabs/dem) 2 • In the optically thin approx- 
imation (where r b s ~ f t), setting the absorber size arbitrarily 
small can spread the overlaid lines over all of the data points (with 
equal spacing on the abscissa, as seen for r b s ^ 1)- This, how- 
ever, is not justified for r J 0.3 and using r D b s = /(l — e ~ T ) 
introduces the compression at large optical depths since r bs — ► / 
asr-> oo [f| Increasing the size of the absorber (e.g. to d a b s ~ 500 
pc) accounts for the four points with d em ^ 2 kpc, but leaves the 



for a uniform source, based upon the limit 



, from T b s to unit y tor a uniform so ' 
imposed by Eq.l2l lo'Dea et all 19941) 



^ In fact, from 

/ > fobs i 

4 Conversely, it is common practice to use Nm, where known from obser- 
vations of the Lyman-a line, in conju nction with f T b B dv, from the 21-cm 
line, to obtain T sp ; n / / in DLAs (See Cumin 2012 and references therein). 

5 Although r b s ~ 0.3 for the vast majority of absorbers. 
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Figur e 2. The peak observed optical depth for the data of IPihlstrom et alj 
120031) [cf. Fig.[T] middle], but with the overlain lines showing how, for a 
given absorber size, the optical depth varies for / = (rf a bs/^om) 2 - Here 
the absorption cross-section is set to d a b s = 100 pc. The values of t, rang- 
ing from 0.001 to 100, are given by the (approximate, for r ^ 1) maximum 
values on the ordinate. 
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Figure 3. As Fig. [2] but for / = d^ s /d em and d a b s = 33 pc (the smallest 
of the bulk values of d cm , i.e. excluding the outlier at d cul fS 10 pc). Again, 
the t — 10 and 100 lines are coincident. 



with the extent of the radio emission being determined by the lobe 
separation. We now explore this, in conjunction with the possibil- 
ity that the large range of emitter sizes may be due to projection 
effects, by extending the sample to also include the non-compact 
sources. 



2.2 Extending to the general population of radio sources 

In discussing the e ffects of radio source size on the strength of the 
21 -cm absorption. ICurran & Whitind j201oh noted that the addi- 
tion of non-compact/unclassified sources did not dilute the anti- 
correlation. Furthermore, from optical, near and far-infrared ob- 
servations, no difference is found b etwe en compact and e xtended 
radio sources jde Vries et al. Hl998l and iFanti et al respec- 
tively). Including th e non-compact/un classified objects for which 
rf em is available (e.g. lGupta et alj2006l) . the Kendall's T two-sided 
probability of the correlation between observed optical depth and 
projected source size arising by chance is P(T) = 0.0044. This is 
in comparison to the P(T) < 0.01 "column density" correlation 
dPihlstrom et alj|2003l) and corresponds to a significance of 2.84<r, 
assu ming Gaussian stat istics (Fig.[4]l[f] 

iGupta et alj feOOfj) state that the incidence of H I absorption 
is much higher for compact sources than for extended sources. 
The histogram in Fig. [4] shows the relative number of detections 
(hatched) and non-detections (unfilled) in each log 10 (d cm ) bin. 
One thing immediately obvious from this is that, while the de- 
tection rates are < 30% in all other bins, between 0.1 and 1 kpc 
the rate exceeds 50%. A peak in the observed optical depth is also 
apparent in this bin (the top panel of Fig. |4). That is, there may 
exist a "resonance" between the absorber and radio source size. 
This could arise if when d a b s <S rf e m the covering factor is low, 
whereas when rf a b s 3> d cm , the chance of a misalignment between 
the absorber and emitter along our sight-line is greater for a smaller 
emission region, indicative of a clumpy absorbing medium. Inter- 
estingly, ~ 100 pc is the scale of the dense st H I clouds found from 
low redshift emission studies dBraunll2~012|) . 

We find the 21 -cm detection rate to be 39% (30/77) at d cm < 1 
kpc, in comparison to 23% (14/60) at d cm > 1 kpc, confirming a 
slightly higher detection rate for the compact sources. However, 
this is elevated by the 0.1-1 kpc bin. Furthermore, if the linear 
extent of the source is correlated with its luminosity, in conjunc- 
tion with t he fact that the rad io and ultra-violet luminosities are 
correlated dAllison et al l2012h . we would expect a lower detection 
rate for the large/more luminous sources, due to ionisation of the 
gas by the active galactic nucleus jCurran & Whitind20T2l) . When 
sources above the critical luminosity (Luv ~ 10 23 W Hz -1 ) are 
removed, 21-cm detection rates in compact obj ects are not signifi- 
cantly higher than for the rest of the population JCurran & Whitind 
l201Ch . 



smaller sources unaccounted for, although, as stated above, there is 
no reason to expect a common absorber size. 

The data are better fit by a covering factor which depends 
linearly on the source size, i.e. / = d a b a /dem, which can trace 
most of the points for a single absorber size (Fig. [3). This sug- 
gests that the geometry of the emission could be dominant along 
a single axis. This may be expected since high resolution imaging 
reveals the sources of emission to usually be two main lobes (e.g. 
iTzioumis et alj2 002: Pihlstrom et al J2003I and references therein), 



2.3 Absorption of a double lobed radio source 

In order to test whether the wide range in observed source sizes 
could, at least in part, be caused by projection effects, we model 
the source size, d em , through the inclination, i, of a small range of 



6 This significance falls t o 1.96a when the no n-detections are included 
(via the survival a nalysis o flLavallevet all 19921) . since these are subject to 
additional effects ICurran & Whitinsdl2010l) " 
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Figure 5. Schematic of the model showing projected linear size, d C m in relation to the lobe inclination, i, and thickness, t owac 320121) . 
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Figure 4. The peak observed optical depth versus the projected linear size 
(where available) for all of the redshifted radio sources searched in 21- 
cm ab sorption to date (compiled in ICurran & Wh iting 2010; Allison et al] 
l2012h .The filled markers/hatched histogram represent the detections and 
the unfilled markers/histogram show the 3cr upper limits to the optical 
depth for the non-d etections at a spectral resolution of 167 kms -1 (see 
ICurran et al]|2013ij) . The circles designate compact objects and the stars 
non-compact/unclassified, with the triangle s showing the high frequency 
peaker galaxies, which, although compact, lOrienti et al.l 1 120061) argue do 
not follow the trend. The line shows the least-squares-fit to the detections. 



intrinsic sizesQThat is, <f om = t+ dcosi, where t is the diameter 
of the lobe and d the deprojected space between the lobe centres 



7 It is the general consensus that, at least some (CSOs & CSSs), compac t 
objects are young and therefore intrinsically small (e.g. iFanti et alii 19951) . 
However, since we are including all of the redshifted associated systems in 



(Fig.H). As described bv lGlowackJ ll2012h . the cross-sectional area 
of the lobes intercepted by a spherical absorber, i.e. the covering 
factor, was recorded when incrementing the inclination. Although 
there may be a preferred absorber size (Sect. 12. 2] , we do not know 
the underlying distribution of the intrinsic properties (such as op- 
tical depth and deprojected size) and so cannot perform a reliable 
statistical fit to these data. Therefore, for various intrinsic optical 
depths, r, the covering factor was converted to the observed optical 
depth via r b s = /(l — e_r ), with the best fit to the trend of the 
detections judged by eye. Different models were tested through the 
variation of several parameters: 

- Varying the relative sizes of the model components: Specifi- 
cally, the size of the absorber in relation to the lobe size, lobe thick- 
ness and spacing between the lobes, as well variation of the relative 
sizes between these three parameters. Generally, an absorber size 
close to or larger than the deprojected spacing between the lobes 
gave the better fit. 

- Varying the lobe morphologies: Several morphologies were 
tested, with the best results being obtained by an ellipsoicjfl 
modelling the domina nt jet morphology of an FRl source 
dFanaroff & Rilevll 19741) . and a "lollipop"|j modelling the faint jet 
plus hot-spot morphology of an FRlI source. 

- Introducing a non-uniform density distribution: In order to 
simulate the clumpy nature of the absorber, the points forming 
the density distribution were generated randomly, where if a value 
(RA n ,<5 n ) is at a radius smaller than the next generated value 
(RAn+i,<5n+i), then absorption occurs at (RA„,<5 n ). As seen from 
Fig. [6] this yields a clumpy distribution in which the density of ab- 
sorbing points falls with radius. Using the non-uniform absorber 
has the effect of smoothing the dependence of the observed optical 
depth on inclination, generally providing a better match to the data 
than a uniformly dense absorber. 

InFigs.[7]and[8] we show distributions obtained for the FRl and FRlI 
approximations, respectively, intercepted by a non-uniformly dense 
absorber of total extent 1 kpc (Fig. [6), which yield ed the r b B — d em 
slopes that best traced the observed distribution (Glowacki 2012). 
From these, just as for the intrinsic optical depths (Sect. I2.lt . the 



which 21 -cm has been detected, it is of interest to see if the inclination of a 
single intrinsic size can account for the observed range in projected sizes. 

8 Which was projected as ellipse that shortened with inclination until form- 
ing a circle at i = 90° . 

9 In which the "stick" shortened with inclination until forming a circle at 
i = 90°. 
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Figure 6. A single realisation of the average coverage with radius for a 
1 kpc non-uniform absorber. The inset shows the density distribution in the 
RA-(5 plane. 
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Figure 7. The observed optical depth versus the projected size (where avail- 
able) for all of the 21-cm detections (see Fig. |4). The overlain lines show 
the variation for a non-uniformly dense absorber of d a b s = 1 kpc, occult- 
ing the elliptical lobe model (FRl), with t = d/10 for d = 0.66, 6.6 and 
66 kpc (the largest of the detections). The maximum value on the abscissa 
shows the deprojected size and on the ordinate the same values of T b s as 
shown in Fig. [2] The inset illustrates the un-inclined 6.6 kpc model in the 
projected sky plane. 



projection of a single intrinsic source size cannot acco unt for the 
whole distribution (as found for the compact objects bv lFanti et alj 
ll99Cb . although a range of sizes in conjunction with a range of 
optical depths does exhibit a similar slope to the r Q b s — d em anti- 
correlation. At large linear extents, there are two non-compact ob- 
jects which are located away from the model "fits", although these 
did not form part of the original sample which led to the perceived 
anti-correla tion (cf. Figs . [2] and [3} . 

Lastlv. |Pihlstromet~afl [2003) model the absorbing gas as both 
a spherical cloud and a disk, centered on the nucleus (the origin of 
the axes in Fig.f5](. In the former case the model results will be the 
same as presented above, since, for a given absorber, the covering 
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Figure 8 

diameter 
d/60. 



As Fig. [7] but for the lollipop model (FRlI), again with a lobe 
of t = d/10, for d = 0.66, 6.6 and 66 kpc and a "stick" width 



factor is solely dependent on the projected area of the lobes inter- 
cepted by the absorber. In the latter case, where the absorption is 
occuring in the pc-scale torus or the kpc-scale galactic disk, with 
a fixed orientation perpendic ular to the jet ax i s, we may expect the 
power law drop as seen by IPihlstrom et al.l d2003l) . although this 
would be with respect to the observed optical depth and not the 
column density (cf. Figs. [T] and |2). However, from the incidence 
of 21-cm detection s in type-1 and type-2 active galactic nuclei, 
ICurran et al.1 J2008I) conclude that the bulk absorption occurs in the 
galactic disk, which must be randomly orientated with respect to 
the torus. This rules out applying a generic model of the jet inclina- 
tion to a disk distribution, although we would generally expect the 
observed optical depth to decrease with decreasing inclination due 
to a reduction in the covering factor (Sect. |2~TV 



3 DISCUSSION AND SUMMARY 

Knowledge of the density of neutral gas in the distant Universe 
is crucial to our understanding of its contribution to the cosmic 
mass density (e.g. lCurranll2010l and re ferences therein) and its re - 
lation to the star formation history (e.g. lHopkins & Beacom|[2006t) . 
Therefore, any dependence of the total column density of the neu- 
tral gas h as profound consequen ces in the interpretation of this 
and, since IPihlstrom et al.l (J2.Q03) noted the column density to be 
anti-correlated with size of the radio source illuminating the H I 
cloud, this relationship has been confirmed several times, although 
no clear physical explanation has ever been provided. 

Upon dissecting the compone nts which make up the estimate 
of the column density, we find, as IPihlstrom et all the correlation 
to be between the observed optical depth and radio source size. The 
former quantity is directly proportional to the covering factor, in the 
optically thin regime, which is itself inversely proportional to the 
source size. Therefore, these quantities are not independent, with an 
observed optical depth-source size anti-correlation expected purely 
from the definition of r b s alone. It is the assumption of a common 
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spin temperature and covering factor for all of the sources, which 
leads to the Nm — d cm anti-correlation. 

Ext ending the sample to i nclude the non-com pact objects, as 
found bv lde Vries etail dl998h ; lFanti et alj ( feOOOh . we find no dis- 
tinct difference between these and the compact sources, with the 
non-compact objects also following the r b s — d cm anti-correlation. 
Also, whi le the 21 -cm detec tion rate may be higher in compact ob- 
jects (e.g. iGupta et"al]l2006l) . this could be dominated by a "reso- 
nance" at d cm pa 0.1 — 1 kpc, in addition to the fact that the larger 
sources may have higher ultra-viol et luminosities, decreasi ng the 
likelihood of detecting neutral gas dCurran & Whitindlioiob . This 
resonance, where the 21 -cm detection rate is significantly higher 
than for other values of d cm , may indicate that the typical cross 
section of cold absorbing gas is of the order 100 — 1000 pc. 

Finally, through a model of a double lobed radio source, where 
the projected linear size, d em , depends upon the inclination, we find 
that the anti-correlation can be best accounted for models which 
approximate either of the two Fanaroff & Riley classifications and 
a non-uniformly dense absorber with a diam eter in the range o f the 
resonant d cm . We also confirm the finding of lFanti et al.l d 19900 that 
the observed range in source extents cannot be accounted for by the 
inclination of a single deprojected size. 
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